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ABSTRACT 

An  analytical  investigation  was  carried  out  to  determine  the 
design  conditions,   for  a  ducted  propeller  unit  operating  in  water, 
providing  optimum  suppression  of  cavitation  when  the  forward  speed, 
depth  and  body  drag  are  given.     The  selection  charts  were  constructed 
in  detail  for  propellers  using  the  NACA  65  series  blades.     The  results 
indicate  the  pumps  of  generally  higher  flow  rates  are  more  desirable 
than  those  that  have  oeen  employed  in  the  few  units  actually  covistructed, 


TABLE  OF  CONTENTS 

Section        Title  Page 
Acknowledgements 
Abstract 

Table  of  Contents 
List  of  Symbols 

I.  INTRODUCTION  1 

II.  DEVELOPMENT  OF  THE  ANALYSIS  4 

III.  COMPUTATIONAL  PROCEDURE  16 

IV.  DISCUSSION  OF  ANALYTICAL  RESULTS  18 

V.  CONCLUDING  REMARKS  23 
References  24 
Figures  25 


LIST  OF  SYMBOLS 
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At        cross-sectional  area  at  entrance  to  guide  vanes  (ft   ) 

C  d       drag  coefficient  of  propelled  vehicle 

D         drag  (lb) 

F  thrust (lb) 

K  cavitation  number  associated  with  minimum  operating  depth, 

at  the  inception  of  cavitation 

Pc         atmospheric  pressure  (lbs /ft    ) 

Pj         local  static  pressure  (lbs /ft   ) 

Pr         vapor  pressure  of  water  (lbs /ft   ) 

2 

Pr         total  pressure  (lbs /ft   ) 

R  radius  of  propeller  from  axis  of  rotation 

S  pressure  coefficient  across  fan  blades 

S  propelled  vehicle  cross-3ectional  area  (ft   ) 

W  velocity  relative  to  rotating  blade  (ft/sec) 

c  chord  length 

d  minimum  operating  depth  (ft   salt  water) 

"t  tangential  blade  spacing  (ft) 

U.  axial  velocity,   absolute  (ft/sec) 

r&  tangential  velocity,    absolute  (ft/sec) 

Z.  number  of  blades 

o<  angle  of  attack  between  flow  direction  and  fan  blade  chord 

/2  flow  angle  between  flow  direction  and  cascade  axis 

^  fan  blade  stagger  angle  between  cascade  axis  and  blade  chord 

T\p  propulsive  efficiency 

©  turning  angle,   the  difference  between  flow  angle  at  entrance 

to  fan  blade  and  exit  from  fan  blade 


X  drag-area  coefficient,   propelled  vehicle  cross-sectional  area 

divided  by  duct  cross -sectional  area,   time3  the  drag  coefficient 

TTd        diffusion  coefficient,    absolute  axial  velocity  at  inlet  to  guide 
vanes  divided  by  free  stream  velocity 

3 

n  density  of  sea  water  (lbs /ft   ) 

CT         solidity,   chord  of  blades  divided  by  tangential  spacing 

(p         flow  coefficient,    rotating  velocity  divided  by  free  stream 
velocity 

~y         turning  coefficient,   tangential  turning  velocity  divided  by  free 
stream  velocity 

u»  angular  rotating  velocity 

Subscripts 

0  free  stream 

1  upstream  of  guide  vanes 

Z.         upstream  of  rotating  fan  blades 
3  downstream  of  rotating  fan  blades 

4-  downstream  of  nozsle 


I.     INTRODUCTION 

The  ducted  propeller  is  of  interest  in  hydrodynamic  applications 
largely  because  it  offers  the  possibilities  of  either  high  loading  at  low 
speeds  or  cavitation  suppression  at  higher  speeds.     The  present  work 
is  concerned  primarily  with  utility  of  the  ducted  propeller  to  suppress 
cavitation  on  the  propeller  blades  as  a  means  of  achieving  reasonably 
efficient  and  quiet  operation  at  high  speeds.     The  several  actual  attempts 
made  to  achieve  this  end  have  not  resulted  in  narked  success  largely 
because  the  principle  wa3  applied  to  problems  where  the  speed,   depth 
or  propeller  design  were  inappropriate.     This  difficulty  is  under- 
standable because,   at  moderate  forward  speeds,   the  selection  of 
propeller  blading,   rotative  speed  and  diffuser  design  must  be  made  in 
a  rather  particular  manner  to  achieve  the  potential  advantage  of  the 
ducted  propeller.     Furthermore  there  existed  no  method  for  selection 
of  the  appropriate  configuration  to  satisfy  given  requirements.     It  is  the 
purpose  of  the  present  analysis  to  develop  a  method  for  finding  optimum 
configurations  for  suppression  of  blade  cavitation  and  to  carry  out 
design  charts  in  detail  for  a  particular  class  of  propeller  blading. 

The  conceptual  model  employed  is  an  axial  flow  ducted  fan  with 
guide  vanes  ahead  of  the  rotating  fan  blades.     Possible  physical 
arrangements  are  shown  in  figure  la,    b  and  c.     The  flow  entering  the 
hydrojet  of  figure  ia  and  lb  is  that  flow  from  the  wake  of  the  vehicle 
being  propelled.     The  flow  entering  the  hydrojet  of  figure  lc  is  the 
free  stream  flow.     This  investigation  was  limited  to  considerations  of 
the  free  stream  flow  entering  the  hydrojet.     Cavitation  was  assumed 
to  occur  on  the  fan  blade;  no  consideration  was  given  to  the  possible 
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occurrence  of  cavitation  due  to  the  geometry  of  the  diffuser  and  duct  or 
of  the  body  being  propelled.     Although  it  is  known  that  cavitation  may 
occur  first  at  the  blade  tip  or  hub,  high  noise  output,   lo~s  of  efficiency 
and  structural  damage  are  largely  associated  with  cavitation  of  the 
blade  surface. 

In  the  development  of  the  analysis,   the  free  stream  flow  velocity 
Uo    was  considered  to  be  ideally  diffused  to  a  lower  velocity    u,   at  the 
entrance  to  the  guide  vanes.     The  flow  through  the  guide  vanes  and 
rotating  propeller  blades  v/a3  analysed  using  airfoil  cascade  theory. 
The  flow  was  treated  at  a  typical  radius   R    from  the  rotating  axis  and 
is  thus  approximately  applicable  to  any  section  of  the  propeller  blade. 
The  guide  vanes  were  considered  to  impart  a  tangential  velocity,    /v-    , 
to  the  flow  entering  the  rotating  blades,   the  rotating  blades  were 
considered  to  turn  the  flow  back  to  an  axial  throughflow,   and  to 
increase  the  total  pressure  by     APt  •     It  was  assumed  that  the  total 
pressure  rise  was  constant  from  hub  to  tip,   that  the  axial  flow  velocity 
was  constant,   and  that  the  radial  velocity  was  negligible.     The  flow 
downstream  of  the  rotating  fan  blades  was  then  discharged  through  a 
nozzle  to  an  exit  velocity    u,^   and  an  exit  pressure     p^    ,    equal  to  the 
free  stream  pressure     po  .     Cavitation  occurred  on  the  suction  side 
of  the  rotating  blade  section  when  the  local  static  pressure   p^  dropped 
to  a  certain  predetermined  value    pv  . 

The  thrust  of  the  unit  was  determined  by  elementary  application 
of  the  momentum  theorem  to  the  flow  upstream  and  downstream  of  the 
ducted  propeller.     Since  steady  motion  was  assumed,   the  thrust  must 
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equal  the  body  drag  and  hence  propeller  requirements  were  related  to 
the  drag  characteristics  of  the  body. 

Performance  data  for  the  rotor  blades  was  extracted  from 
references  (1)  and  (2)  for  the  NAC.A  65  series  airfoil  in  cascade.     These 
results  were  employed  to  obtain  a  typical  3et  of  selection  charts, 
applicable  only  to  these  blades,   to  illustrate  use  of  the  analysis  to  select 
appropriate  blade  sections  for  a  given  forward  speed,   depth  of  operation 
and  body  to  be  propelled. 


II.  DEVELOPMENT  OF  THE  ANALYSIS 

The  thrust  of  a  ducted  propeller  is  derived  from  the  change  in 
momentum  of  the  fluid  flow  thru  the  unit.     Thi3  change  in  momentum 
ia  accomplished  by  the  transfer  of  energy  between  the  fluid  and  the 
rotating  propeller  blades  resulting  in  a  total  pressure  increase  which 
is  converted  to  kinetic  energy  in  the  nozzle.     The  mass  flow  thru  the 
duct  is  described  by     o  u,  A| ,   where    o    is  the  density  of  sea  water,  u.( 
is  the  absolute  axial  flow  velocity  upstream  of  the  guide  vanes,   and  A| 
is  the  cross -sectional  area  of  fluid  flow  at  the  inlet  to  the  guide  vanes. 
By  continuity  for  fluid  flow  and  since      o    is  constant  for  an  incom- 
pressible fluid,   the  mass  flow  upstream  and  downstream  of  the  hydro- 
jet  is  also     c  u,  Ai  •     The  change  in  momentum  thru  the  ducted  pro- 
peller is  therefore  the  mass  flow  times  the  change  in  velocity  between 
upstream  and  downstream  states,   and  is  equal  to  the  thrust  of  the  hydro 
jet 

p  =   <=u,A((uh  -Uo)  (2.1) 

This  thrust  must  equal  the  drag  of  the  unit  being  propelled  in 
order  to  maintain  a  forward  speed  of    u0  .     Drag  is  expressed  as 

D  =  ±  f  Uo'CpS'  (2.2) 

where     Co     i3  the  drag  coefficient  defined  so  as  to  include,   in  addition 
to  usual  aerodynamic  drag,   the  pressure  forces  over  the  streamlines 
separating  the  flow  passing  through  the  pump  from  that  passing  around 
it.      S        is  the  cross -sectional  area  of  the  body  being  propelled.     From 
equations  2.  1  and  2.  2  we  have 


_^_ 


'Ao 


=  i(^')i^7)+l  «2-3> 


Propulsion  efficiency  is  defined  as  the  thrust  work  divided  by 
the  work  input  and  is  expressed  as 


which  reduces  to 


np=H¥To 


From  equations  2.  3  and  2.  5  it  is  seen  that  two  major  dimen- 
sionless  parameters  appear  that  are  important  to  efficiency  and  to 
the  size  requirements  of  the  hydrojet  relative  to  the  size  of  the  vehicle, 
These  are  defined  as 


A=  4s-  •      7Td  =   -£-  (2-6) 


where     A     is  a  combined  drag  and  cross-sectional  area  coefficient 
and   7Td  is  a  diffusion  coefficient  from  free  stream  to  the  guide  vanes. 
The  interrelations  of      A      »     7Td    and    t\f>     are  shown  in  figure  3. 

To  analyze  the  guide  vanes  and  rotating  fan,   the  airfoil  cascade 
theory  was  applied.     An  imaginary  cylindrical  cro3s-section  of  the 
blades  is  taken  with  axis    coinciding    with  the  axis  of  rotation.     The 
imaginary  cylinder  is  assumed  to  have  radial  extent    dR     and  a  radius  R 
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and  then  ia  imagined  to  be  unrolled.     A  row  of  parallel  blade  cross  - 
sectional  cuts  results  which  extends  on  to  infinity  if  the  pattern  is  re- 
peated indefinitely.     The  flow  problem  then  is  reduced  from  a  three- 
dimensional  to  a  two-dimensional  problem.     The  rotational  velocity  of 
the  rotor  blades  is  represented  by  a  tangential  velocity  ujR    perpendic- 
ular to  the  cascade  axis.     An  observer  can  assume  that  he  is  fixed  in 
a  coordinate  system  relative  to  the  rotating  blade  sections  or  in  an 
absolute  system  with  respect  to  the  guide  vanes.     The  difference  in 
the  two  coordinate  system  then  is  only  in  the  rotating  velocity,    ujj^  . 

The  relative  spacing  of  blade  sections  is  described  by  the 
solidity,    o~  ,   and 

cT  -     c/t  <2-7> 

where     C     is  the  blade  section  chord  length  and     t     is  the  tangential 
spacing  equivalent  to     2  7T  R    divided  by  the  numoer  of  blades.     The 
solidity  parameter  is  one  of  the  important  design  quantities  which 
effects  the  performance  of  the  airfoil  sections  and  the  occurrence  of 
cavitation. 

In  figure  Z  is  shown  a  typical  row  of  guide  vanes  and  rotor  blade 
sections  obtained  by  the  unrolling  of  a  cylindrical  cross -section. 
Reference  stations  are  designated  as  (0)  to  indicated  free  stream  flow, 
(1)   to  indicate  the  flow  entering  the  guide  vanes,    (2)  as  exit  from  the 
guide  vanes  and  entrance  to  the  rotating  fan,    (3)  as  exit  from  the 
rotating  fan,   and  (4)  as  the  exit  from  the  nozzle.     The  angle  between 
the  blade  section  chord  and  cascade  axis  i3  designated  the  stagger 
angle,  J    .     The  direction  of  flow  relative  to  the  rotating  blade  sections 
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is  established  by  the  inlet  flow  angle,    Sz  ,   measured  from  the  cascade 
axis  to  the  direction  of  flow.     Similarly,   the  exit  flow  angle    /S 3    ,   is 
measured  from  the  cascade  axis  to  the  direction  of  flow  relative  to  and 
leaving  the  rotating  blades.     The  difference  between  the  inlet  flow  angle 
and  the  exit  flow  angle  is  the  turning  angle,    0    ,   thru  which  the  flow  is 
turned  by  the  rotating  blades. 

The  flow  through  the  ducted  propeller,   absolute  and  relative  to 
the  rotating  blades,   can  be  represented  by  a  velocity  triangle  as  shown 
in  figure  2.     The  axial  flow  velocity  was  assumed  constant  therefore 
U,     =    U3      .     The  flow  velocity  relative  to  the  rotating  fans  is  obtained 
by  subtracting  the  rotating  velocity,    oo  R    ,   from  the  absolute  velocity, 
Ux   ,   entering  the  blades  sections,    giving  the  inlet  flow  velocity,    W^   , 
relative  to  the  rotating  blades.     The  relative  flow  velocity  leaving  the 
fan  is    W3     ,   and  is  converted  to  an  absolute  velocity,    u.^    ,   by  adding 
the  rotating  velocity,  u;  R  ,   to  the  relative  velocity.     The  tangential 
velocity,  AJ-  ,    imparted  to  the  fluid  by  the  rotating  blades  was  assumed 
to  be  equal  to  the  tangential  velocity  imparted  by  the  guide  vanes, 
resulting  in  the  simplified  velocity  triangles  as  given  in  figure  2. 

To  describe  the  pressure  distribution  across  the  rotating  blade 
section,   the  pressure  coefficient,    S    ,   is  used  as  given  in  references 
( i)  and  (2)  where 

c  =    P<T  ~*  P)l  t?   8) 

b     tjvtf  (2*8) 

the  total  pressure  entering  the  rotating  blade  section,   minus  the  local 
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static  pressure  at  any  point  across  the  blade  section,   divided  by  the 
inlet  dynamic  pressure.     When  the  local  pressure,     P£     ,   is  equal  to 
some  predetermined  pressure    Pv    ,   cavitation  is  assumed  to  occur. 
The  value  of    Pv    may  be  taken  as  the  local  vapor  pressure  of  the 
water  or  a  more  suitable  figure  if  knowledge  of  local  conditions 
permits  it.     Thus  for  given  inlet  conditions  there  will  be  a  maximum 
value  of  the  pressure  coefficient,     S  ha\  »   for  operation  without 
cavitation.     To  relate  the  free  stream  conditions  to     Sn^x      »    a 
cavitation  number    K     ,   is  employed  where 

.      Po-Pv 


ft   ^    1st L*_  (2.9) 


the  free  stream  static  pressure  minus  the  cavitation  pressure  of  sea 
water,   divided  by  the  free  stream  dynamic  pressure.     It  will  be  shown 
that  the  value  of     ft     determines  the  minimum  operating  depth       6 
corresponding  to    P0     for  a  given  operating  velocity,     U0    .     Part  of  the 
problem  of  this  analysis  is  to  relate  the  cavitation  number     K   ,   to  the 
pressure  coefficient,      S  mwc     »    where      5  rift*,     is  known  from  experi- 
mental data. 

Two  other  dimensionless  parameter s  employed  for  this  analysis 
are  the  flow  coefficient     <p     ,   and  turning  coefficient  ^    where 


9  -  ~^  <2.  10) 


tv 


Y--~  (2.1D 
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thus  relating  the  rotational  velocity  and  the  turning  velocity  to  the  free 
stream  conditions. 

By  the  assumption  that  no  losses  occur,   the  total  pressure  rise  , 
APT»    across  the  propeller  can  be  obtained  by  subtracting  the  free  stream 
total  pressure  from  the  exit  total  pressure  so  that 

APT    -     Pmt    -   Pot  (2.  12) 

The  total  pressure  is  equal  to  the  sum  of  the  static  and  dynamic 
pressures  so  that 

PHT  -   PH    +    ^y  U*  (2.  13) 

and 

Pot    -   Po     +    -k   f**  (2.  14) 

Since  the  static  pressures  at  entrance  and  exit  are  assumed  to 
be  equal, 

AP-  i  ?  ^  [  af  -  0  <*•  "I 

However,    since  no  losses  occur  and  the  axial  flow  velocity  was 
assumed  constant,   we  obtain 

PHT-  P,T  =  P5    +    i   ?  ^  (2.16) 


Pot  =  Pit    -  R    +    -i  J  U.,*  (2.  17) 


? 
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and 

APT  =    P,   -  P,  (2.  18) 

Thus,    Bernoulli's  equation  relative  to  the  rotating  blades  can  be  applied 
to  obtain  the  total  pressure  rise  across  the  hydrojet.     Therefore, 
relative  to  the  rotating  blades 

K  +"   i  f  W*    =     P3  +  i  f  W3*  (2.1)) 

From  the  velocity  triangle  of  figure  2,   it  is  seen  that  the  relative 
velocities  are 

W2*   -  U*    4-     (at  -+  ujR)  (2.20) 

W*    =  ^    +     HR)  (2.21) 

where  again  the  axial  velocities     u.,      and      u^      are  equal.     Relative  to 
the  guide  vane,    Bernoulli's  equation  would  be 

P,  +  t  ?  u«t    =    P2   +  t  f  as-  (2.  22) 


where 


U*  =     Ah*  +    a,  (2.  23) 


By  combining  equations  2.  18  thru  2.  23  Euler's  turbine  equation  is 
obtained  as  follows: 

AP  =   P3   ~P,    =     ?  ^R  rtf  (2.24) 
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Therefore  from  equations  2.  15  and  2.  24  we  obtain 

A?T  =  ?  uu  R  rJ-   =    I  f  ^  [  (~;)    -  j)  (2.  25) 

or  in  terms  of  the  flow  coefficient     <p    and  turning  coefficient      ^ 

^T=(^f-I  (2.26) 

When  thrust  and  drag  were  equated  it  was  found,   equation  2.  3,   that 

ig-ik  +  l  U-27) 

Therefore 

Thi3  important  relationship  matches  the  operating  characteristics  of 
the  ducted  propeller,     dp    and   "li/-     ,    with  the  drag  of  the  vehicle  being 
propelled,   the  relative  size  of  the  vehicle  and  duct  cross-section,    and 
the  performance  of  the  entrance  diffusor.     In  a  sense,   the  product 
is  a  measure  of  the  power  requirements  of  the  ducted  propeller. 

Associated  with  the  rotational  speed  and  turning  velocity  is 
the  turning  angle,    0    ,   which  is  the  difference  in  the  flow  directions  at 
the  entrance  and  exit  relative  to  the  rotating  blades,   therefore 

tan    9        i    .    +.      /o  t""j  (2.29) 

But  from  the  velocity  triangle  of  figure  2  it  is  seen  that 
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-to-n/3^-        ^  (2.30) 


tan/5j=    -^-  (2.31) 


By  substitution  and  by  dividing  numerator  and  denominator  of  equation 
2.  29  by  the  square  of  the  free  stream  velocity,   the  turning  angle  can 
be  determined  to  be 

tM  e  "■   ycAyj  +  TV  (2-  32> 

From  equation  2.  32  a  value  of  &  can  be  calculated  which  then  must 
be  matched  with  experimental  performance  data  for  the  particular 
type  airfoils  being  employed  for  the  propeller  olades.     Associated  with 
this  turning  angle  will  be  a  pressure  coefficient,     S      ,   which  can  be 
allowed  to  have  a  maximum  value  when  the  local  pressure  across  the 
blades  is  equal  to  the  pressure     pv    . 

from  the  definition  of  the  pressure  coefficient  we  have  that 


But 

h  =   Po  +  i   ?  [^  -   ("N-  ^.)   3  (2-  34) 


therefore  by  dividing  numerator  and  denon  i    i  tor  .  y      |  o  u3 
we  obtain 


L 
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5-1+     ^^-^-^ ,    ,    ,    _, (2.35) 


When  the  local  pressure  equals  the  pressure  Pv  we  have   SnAX   and  since 
the  cavitation  number   K    was  defined  as 

we  can  determine  the  cavitation  number  as  a  function  of  the  ducted 
propeller  performance  and  operating  characteristics 

k--  is^-oBf+Yf  +  Ttf]  +[n;+r-i]      (*•»*> 

From  the  definition  of  the  cavitation  number  it  is  seen  that  for  a  given 
operating  velocity    U&    the  minimum  operating  depth   d     is  determined 
since  the  free  stream  static  pressure     R>    is  equivalent  to  the  pressure 
at  the  depth    d     plus  the  atmospheric  pressure.     The  expression  relating 
^    ,    u^o    and    d     ,   assuming  the  vapor  pressure  of  water  is  zero,    is 


d«    &£  ft    -    JL-  (2.37) 


and  is  presented  as  a  nomograph  in  figure  19.     Therefore,   from  equations 
2.  36  and  2.  37  the  operating  depth  and  velocity  can  be  determined  as  a 
function  of  the  performance  of  the  rotating  blades,   the  diffusion  and 
power  requirements  for  the  ducted  propeller.     To  operate  at  shallow 
depths  with  high  velocity  it  is  seen  that  the  pressure  coefficient  should 
have  a  minimum  value  for  the  maximum  suppression  of  cavitation. 


-14- 

This  in  turn  requires  that  the  pressure  coefficient  should  have  as 
small  a  value  possible  which  is  equivalent  to  saying  that  the  velocities 
across  the  blade  sections  must  be  minimized.     Therefore,   to  suppress 
cavitation  there  will  be  optimum  values  for  all  of  the  dimensionless 
parameters  discussed  which  will  match  with  the  actual  performance  of 
the  blade  sections  as  given  in  references  (1)  and  (2)  for  airfoils  in 
cascade. 

In  this  section,   the  important  relationships  have  been  derived. 
In  the  next  section  the  computational  procedure  carried  out  to  solve 
the  relationships  and  to  match  the  solutions  with  the  experimental 
performance  data  for  the  65  series  airfoil  sections  in  cascade  will  be 
given.     A  summary  of  the  important  relations  are  as  follows: 


rip*  -j— - (2.5) 

A     f   | 


to,n  6  =       r       "^ —  (2.32) 


K-(SMft*-|)[((p-R}^-h]T;]  +  [V-hYX'-|]  (2.36) 


15 


ft 


(Pv  =  0) 


(2.37) 
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III.     COMPUTATIONAL  PROCEDURE 

From  equations  2.5  and  2.  28  it  is  seen  that  for  fixed  diffusion 
coefficient    7Tj    and  drag-area  coefficient    \  ,   the  propulsive  efficiency 
flp  ,   and  the  product  of  the  flow  coefficient     <p     and  turning  coefficient"^ 
are  fixed.     Then  using  equation  2.  32,   the  turning  angle     G  can  be 
determined  for  a  range  of  values  of    (D    for  a  series  of  constant  values 
of    TXa    anci     A   •     This  was  done  for  values  of    1S±     between  .  4  and 
1.  1    and  for  values  of   }\   between  .25  and  .75.     It  was  necessary  then 
to  match  the  computed  values  of  the  turning  angle  with  the  performance 
data  for  65  series  airfoils  in  cascade  in  order  to  determine  the 
pressure  coefficient,  SH/rt. 

The  pressure  coefficient  5  max,  is  a  function  of  the  solidity,    CT     , 
angle  of  attack,    oc.    ,    stagger  angle,  J    ,   camber  and  thickness  dis- 
tribution.    The  data  presented  in  reference  (1)  gives    5     and     &      as  a 
function  of     o<     for  values  of  solidity  equal  to  .5,    .75,    1.00,    1.25 
and  1.  50,   and  for  values  of  flow  inlet  angle,  ^    ,   equal  to  30   ,   45    , 
60      and  70    .     The  data  is  also  presented  for  various  airfoils  of  the 
65  series  with  same  thickness  distributions  but  different  values  of 
camber.     In  reference  (2),   data  is  presented  for  airfoil  sections  with 
the  camber  reversed  in  an  attempt  to  reduce  velocities  across  the 
blade  sections. 

Fixing  the  inlet  flow  angle    f$2    ,   and  varying  the  angle  of  attack 
c<    >   requires  that  the  stagger  angle  be  varied.     However  for  a  fixed 
blade  row  the  stagger  angle  must  remain  fixed  requiring  that  the  inlet 
flow  angle  vary  with  the  angle  of  attack.     Therefore  to  use  the  data 


presented  in  references  (i)  and  (2)  it  was  necessary  to  cross  plot  the 
data  to  obtain  data  for  fixed  airfoil  geometry.     The  cross  plotting  was 
carried  out  for  3tagger  angles  of  45     and  55    .     Figure     5      is  an 
example  of  the  results  for  the  65    4i0  airfoil  in  cascade  with  a  solidity 
of  ..0,    stagger  angle  55    ,   and  with    5n«ar:d    Q    plotted  vs    c<    . 

The  computed  values  of  turning  angle,     0    ,    were  then  entered 
into  the  cross  plots  to  determine  the  maximum  value  of  the  pressure 
coefficient,    5  ^x.    ,   corresponding  to  the  values  of    7^     .    \    ,     <p 
and   ^    .     Having     S^AX    ,   the  cavitation  number    K    was  calculated  and 
plotted.     Figure    10     shows  an  example  of  such  a  plot  for    A        .5, 
o        55    ,    CT      1.0  where  the  curves  represent  values  of  K    and     <p    for 
fixed  values  of    7\j  .     A  cross  plot  was  then  made  for  constant  values 
of  the  cavitation  number    K     as  a  function  of    1X&    and    <D  .     Figure    13    is 
an  example  of  the  contour  plots  of   K     obtained  and  since  a  minimum 
value  of    K      is  required  for  the  maximum  suppression  of  ca    itation 
the  design  point  for  the  given  conditions  is  the  center  of  the  contour 
representing  the  minimum  value  of  the  cavitation  number. 
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IV.     DISCUSSION  OF  ANALYTICAL  RESULTS 

To  determine  the  effect  of  the  solidity,   three  cavitation  number 

contour  plots  were  made  for  values  of    0"    =  l.C,    1.25  and  1.50  where 

the  stagger  angle  was  fixed  at  55     and     A-   =£_S  was  fixed  at  .  5  . 

Ai 

Figures  13,    15  and  16  are  the  result  and  clearly  show  that  solidity 
should  have  a  value  of  approximately  1.  25.     Similar  plots  for  values  of 
(X       equal  to  .75  and  .5  were  not  made  because  by  inspection  of  the   Shax 
and   ©  vs    &<      plots  (not  shown)  the  values  of     S  m*x      are  somewhat 
higher  than  for  the  values  shown  in  figures  5,    6  and  7. 

The  effect  of  stagger  angle,    which  is  another  important  design 

parameter,   can  be  determined  by  inspection  of    S^and   Q     plots.     The 

So 
mkk    for  stagger  angles  less  than  45     are  larger  near  design 

point  (minimum     5ni\n,  )  but  will  give  a  larger  spread  of  values  of     <p 

for  minimum  values  of    K      •     For  stagger  angles  greater  than  55     the 

values  of       Snnx       are  again  larger  and  also  the  spread  of  values  of    <p 

is  more  limited.     In  figures  ii  and  13  a  comparison  for  stagger  angles 

of  45     and  55     can  be  made  with  all  other  conditions  remaining  fixed. 

The  increase  to  stagger  angle,    J    ,   from  45     to  55     has  shifted  the 

design  point  slightly  to  the  right,   decreased  the  minimum  value  of  the 

cavitation  number,   and  increased  the  area  inclosed  by  values  of 

constant     K    •     Any  further  increase  to  the  stagger  angle  can  be 

expected  to  shift  the  design  further  to  the  right  and  increase  the  minimum 

values  for    K    .     The  physical  significance  of  shifting  the  design  point  to 

the  right  is  that  larger  values  of    <p    can  be  employed  which  means  that 

the  rotational  velocity  can  be  increased  for  a  given  /alue  of    M,0   ,   and 
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that  the  loading  of  the  blades  is  decreased  by  a  decrease  in  the  angle 
of  attack.     The  optimum  value  for  the  stagger  angle  appears  to  be 
approximately  55     for  the  airfoils  considered. 

The  effect  of  drag-area  coefficient,  A     ,   is  quite  pronounced 
in  shifting  the  design  point.     Calculations  were  made  and  plotted  for 
values  of    \    -  .25,    .50  and  .75  .     The  value  of    A    =  0.25  is  generally 
appropriate  for  a  good  torpedo  shape.     The  resulting  cavitation  number 
contour  plots  are  shown  in  figures  12,    i3  and  i4.     In  addition,   constant 
values  of    K   are  cross  plotted  to  figure  3  to  show  the  effect  of     A     on  K 
as  well  as  propulsive  efficiency,    H,p    .     From  the  figures  it  is  seen 
that  increasing    A      shifts  the  design  point  to  the  right,   increasing    #    , 
and  up,   increasing  the  design  value  for  the  diffusion  coefficient     Tfj  . 
It  also  radically  increases  the  minimum  values  of     K      that  can  be 
obtained,   thus  restricting  the  operating  velocity  and  minimum  depth  of 
operation.     Therefore,   for  the  maximum  suppression  of  cavitation, 
the  drag-area  coefficient,  A     ,    should  be  as  small  as  possible.     This 
requires  a  small  value  at  the  drag  coefficient,    Cd     ,    and  a  small  value 
for  the  ratio  of  cross-sectional  area  of  unit  being  propelled  to  the  cross- 
section  area  of  the  duct.     If  the  drag  at  the  unit  is  high,    a  coxnparatively 
large  ducted  propeller  will  be  required.     The  value  of      A     therefore 
has  a  marked  effect  in  determining  the  maximum  size  of  a  vehicle  for 
propulsion  by  ducted  propeller. 

To  reduce  the  problem  of  obtaining  diffusion  ratios  other  than 
1.  0,   a  maximum  value  of     A      is  desirable  since  the  design  value  of  7T<j 
increases  with     A     .     In  order  to  obtain  a  larger  range  of  values  for 
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the  flow  coefficient    <p   figures   12  through  14  indicate  that  a  minimum 
value  for    A    is  required.     Since    ip  <^>  ft     for  given  speed  and  angular 
velocity,  large  range  in     ^     is  desirable  to  permit  a  large  hub  ratio 
and  hence  a  small  outer  diameter.     From  figure  3  it  is  evident  that 
decreasing    A     will  increase  the  propulsive  efficiency    Y\p  .     In  view 
of  the  above  considerations,   a  low  value  for    A      is  again  more 
desirable  than  a  high  one. 

The  effect  of  increase  in  camber  of  the  blade  sections  is  to 
increase  the  velocities  across  the  blade  surfaces  and  thus  increase  the 
pressure  coefficient,    Stym.     Calculations  and  plots  were  made  for  the 
65-410  and  65-810  airfoil  sections  for  comparison  where  stagger  angle 
and   A    were  fixed.     The  results  are  shown  in  figures   13  and  17  and 
show  that  increase  in  camber  from  4    /o  to  8     /o  shift3  the  design  to 
the  left,   decreasing  the  design  value  for   <p   .     The  minimum  value  for 
the  cavitation  number,  K    ,   has  not  been  appreciably  affected  but  the 
areas  enclosed  by  constant  values  of    K    have  been  reduced  considerably, 
reducing  the  range  of   (D  values  available  for  a  given  design  value  for 
K    .     The  plots  of    S„A  and     &     vs     o<     for  airfoil  sections  with  greater 
cambers  show  that     SM^increases  considerably  and  that  the  range  of 
useful  values  of     (p    will  be  reduced  further. 

In  figure  18,   the  resulting  contour  plot  for  the  65-  (8AtI8b)lO 
airfoil  sections  is  shown  to  point  out  the  effect  of  reversing  the  camber 
distribution  to  reduce  local  velocity  ratios.     Although  this  decreases 
the  minimum  value  of    K     ,   the  range  of  values  for     <p      at  fixed    7Td     , 
is  decreased  severely.     Thus  it  does  not  appear  advantageous  to  reduce 


local  velocities  over  blade  surface  ijy  this  tech.. 

To  illustrate  how  data  giver,  in  the  figures  can  be  employed  to 
determine  esign  conditions  for  the  maximm  on  of 

cavitation,   assume  that  a  unit  is  to  be  propelled  at  a  velocity  of  75  ft/sec 
(45  ..riots)  and  is  to  operate  at  a  depth  of  ten  feec.  m  the  previous 

discussion  it  wa^  seen  that  a  minimum  value  of  the  drag-area  coefficient 
is  required.     Therefore,   design  considerations  for  the  unit  will  be  for 
as  low  a  drag  coefficient,    Co    ,   as  possible,   and  for  as  favorable  a 
cross-sectional  area  ratio  as  possible  consistant  with  size,    weight  and 
space  available  to  the  unit  being  propelled.     Assume  that  the  Lest  value 
oi'     A    that  can  be  obtained  is  C.5  .     From  figure  i9  it  is  found  that  the 
maximum  allowable  value  for  the  cavitation  number  is  approximately 
0.5  .     For  any  value  greater  than  0.5,   cavitation  can  be  expected  to 
occur.     For     K   =  0.5,   figure  i5  for  the  65-410  i-iade  sectio.:  appears 
to  have  the  largest  area  enclosed  by  a  cavitation  num  .or  oi"  0.5  .     It 
gives  a  large  range  for  the  /alues  of    (p  requires  a  diffusi. 

coefficient,  TTd    i    o£  0.61  .     The  solidity  and  stagger  angle  in  thiu  case 
are  1.25  and  55     respectively.     The  center  of  the  contour  representing 
the  minimum  values  of     K      would  be  designated  as  the  design  point 
for  the  mean  radius     R     between  hub  and  tip.     This  value  for    <p    -  i.  15 
then  would  fix  the  optimum  angular  velocity  uj    and  thus  fixes  the 
values  for    <p     for  changes  in  radius     %.     from  hub  to  tip.     The  ratio 
of  the  maximum  and  minimum  flow  coefficient  contained  within  the 
boundaries  for    K    =0.5  and    TTa  =  0.61  is  approximately  i.  3  .     This 
indicates  that  a  relatively  large  si^e  iiuuj  will  be  necessary,   the  actual 
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size  being  fixed  by  the  value  of  cross -sectional  area  decided  upon  in 
fixing     ,\    . 

From  figure  3,   the  propulsion  efficiency  can  be  deter  mined. 
Then  from  the  propulsive  efficiency  and  drag  coefficient,    the  power 
requirements  can  be  determined  by  replacing  the  thrust  work  by  drag 
and  solving  for  work  input  in  equation  2.  4  .     Thus  a  design  criterion 
has  been  established  for  the  maximum  suppression  of  cavitation. 

If  optimum  diffusion  coefficient  value3  could  not  be  obtained,    if 
greater  velocities  were  required,   or  if  the  value  of  the  drag-area 
coefficient    A    were  unfavorable,   it  would  be  necessary  to  increase  the 
minimum  depth  of  operation  and  hence  raise  the  free  stream  pressure 
P.    . 
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V.     CONCLUDING  REMARKS 

An  attempt  has  been  made  to  determine  what  design  and  operating 
criteria  are  necessary  to  select  a  ducted  propeller  propulsion  unit  with 
the  maximum  suppression  of  cavitation.     The  detailed  selection  charts 
were  limited  to  the  NACA  65 -series  airfoils  for  propeller  blades.     The 
technique,   however,   is  general. 

The  most  sensitive  parameter  appears  to  be  the  drag-area 
coefficient,    \  .     This  quantity  should  have  as  low  a  value  as  possible 
in  order  to  increase  the  range  of  propeller  operation  without  cavitation. 
The  diffusion  coefficient    TVj  is  also  a  significant  parameter  and  it  is 
an  interesting  result  that  an  optimum  retardation  of  the  stream  exists; 
the  stream  may  be  slowed  down  too  much  as  well  as  too  little. 

The  flow  coefficient  parameter,    <p  ,    was  found  to  have  a  rather 
low  optimum  value,   requiring  relatively  3low  rotating  velocities  or 
through  flow  velocities  much  larger  than  conventional.     In  addition  it 
was  found  that  the  range  of  values  for   dp    without  cavitation  will  be  some- 
what limited.     The  design  criteria  for    (p    are  largely  determined  by 
the  blade  sections  being  employed  for  the  propeller,    requiring  low 
velocities  across  the  blades  with  high  loading  of  the  blades. 

For  the  blade  sections  investigated,  it  is  found  that  there  are 
optimum  values  for  the  solidity  (T  ,  the  stagger  angle,  o  ,  and  the 
camber.  The  optimum  values  are  those  that  will  provide  minimum 
values  of  the  maximum  pressure  coefficient,  SMPkA  ,  over  as  large  a 
range  of  angles  of  attack  as  possible,  and  with  limit  on  the  maximum 
values  of  the  turning  angle. 
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a)  Wake  flow  entering  hydrojet 


b)  Wake  flow  entering  hydrojet  driven  by  conventional 
propeller  snaft  . 


c)  Free  stream  flow  entering  hydrojet 


Figure   i.     Schematic  of  possible  hydrojet  propulsion  configurations 
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Figure  2.     Schematic  of  guide  vanes  and  hydrojet  blade 

sections  in  cascade  and  the  associated  velocity 
triangle. 
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Figure  3.        Diffusion  coefficient,     7Td    ,    versus  propulsion 

efficiency    qp  ,   for  constant  drag-area  coefficient,  A    . 
Constant  values  of    K     are  shown  for  the  combination 
<X    =  I.E5       ,    3     =  55°   and   NACA  65-410  blade  section. 
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Figure  10, 


Cavitation  number   K    vs.    flow  coefficient    fy   for 
constant  diffusion  coefficient  TTd  for  the  cascade 
combination      &"    =   1.0,    $    =  55°,    A    =  .5      and  NACA 
65-410  blade  section. 
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1.2 


Figure   11.     Diffusion  coefficient    Hi  vs.    flow  coefficient   <P   for 

constant  values  of  cavitation  number    K    for  the  com' 
bination     A   =  .  5,      CT  =   1.00,     3     =  45°    and   NACA 
65-410  blade  section. 
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Figure   12.     Diffusion  coefficient    Md  vs    flow  coefficient    <p    for 
constant  values  of  cavitation  number    K     for  the  com 
bination    X    =.25,     0*  =  1.  0,  J  =  55°    and  -N  AC  A    J 
65-410  blade  section. 
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Figure   13.     Diffusion  coefficient   7Td  vs    flow  coefficient    <p    for 
constant  values  of  cavitation  number    K     for  the 
combination     \    =  .  5,     C  -  1.00,    3    =  55°  and  NACA 
65-410  blade  section. 
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Figure   14.     Diffusion  coefficient     7\&  vs    flow  coefficient    <p    at 
constant  values  of  cavitation  number    K     for  the 
combination     A    =    .  75,     C  =   1.00,   3    =  55°  and  NACA 
65-410  blade  section. 
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Figure  15.     Diffusion  coefficient  Tm  vs      flow  coefficient     <p     for 
constant  values  of  cavitation  number    K    for  the  com 
bination    A    =  .  5,     0"    =  1.25,    3    =  55°  and  NACA 
65-410  blade  section. 
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Figure   16.     Diffusion  coefficient    7\<j  vs    flow  coefficient    (p     for 
constant  values  of  cavitation  number    K     for  the 
combination     X    =  .  5,     0~    =  1.5,    3    =  55     and  NACA 
65-410  blade  section. 
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Figure   17.     Diffusion  coefficient    7\a  vs    flow  coefficient  (p   for  constant 
values  of  cavitation  number  K    for  the  combination  A  =  .  5, 
<r     =  1.00,    3    =  55°    and   NACA  blade  section  65-810. 
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Figure  18.     Diffusion  coefficient  7\±  vs    flow  coefficient    9  *or 

constant  values  of  cavitation  number  K     for  the  com- 
bination   A     =  .5,    CT    =  1.  00,     3     =55°   and  NACA 
65-  (8A2  l8b)10    blade  section. 
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Figure  19.     Nomograph  relating  the  cavitation  number    K   ,    free 
stream  velocity    U0  ,    and  depth  of  operation    d   . 
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